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Much of the current experimental efforts for detecting Majorana zero modes have been centered on 
probing the boundary of quantum wires with strong spin-orbit coupling. The same type of Majorana 
zero mode can also be realized at crystalline dislocations in 2D superconductors with the nontrivial 
weak topological indices. Unlike at an Abrikosov vortex, at such a dislocation, there are no other low- 
lying midgap states than the Majorana zero mode so that it avoids usual complications encountered 
in experimental detections such as scanning tunneling microscope (STM) measurements. We will 
show that, using the anisotropic dispersion of the tig orbitals of Ti or Ta atoms, such a weak 
topological superconductivity can be realized when the surface 2DEG of SrTiC >3 or KTaOs becomes 
superconducting, which can occur through either intrinsic pairing or proximity to existing s-wave 
superconductors. 


I. INTRODUCTION 

Non-Abelian braiding statistics of well-separated Ma¬ 
jorana zero modes can provide one simpler means for re¬ 
alizing topological quantum computations mm- Partly 
motivated by this, the search for Majorana zero modes 
in nature has become one of the central and challeng¬ 
ing issues in condensed matter physics in last few years 
GSB|. It has been recognized in recent years that super¬ 
conductivity in a system where the spin-orbit coupling 
and the Zeeman field co-exist can be topologically non¬ 
trivial even with conventional s-wave pairings, giving rise 
to Majorana zero modes in topological defects EHS]. Al¬ 
though much of the experimental investigations into this 
physics have focused on ID systems, e.g. Refs [9) [TO] . 
the results are known to hold in two-dimentional (2D) 
systems as well. One promising arena to realize such 2D 
systems is the two-dimensional electronic gases (2DEGs) 
formed at the boundaries of perovskite transition metal 
oxides El- Notable examples of such 2DEGs include 
the (001) interface between SrTiOs and LaA 103 [12], the 
surfaces of SrTiOa (TTij [Dj, and KaT’aO.a pT]- All of 
them possess nonzero Rashba-type of spin-orbit coupling 
because of the lack of inversion symmetry at the bound¬ 
aries. In addition, ferromagnetism is a frequent feature 
of these 2DEGs 117)1 - 1117] even though spin-orbit coupling 
usually tends to suppress spin alignment. Given that in¬ 
trinsic superconductivity has been observed in many of 
these 2DEGs [T5H25] , not to mention the possibility of in¬ 
ducing superconductivity on the surface 2DEGs through 
superconducting proximity effect, one naturally asks the 
question whether topologically protected Majorana zero 
modes can be achieved in the oxide 2DEG. 

One major difficulty in investigating the topology of 
such oxide 2DEG superconductor is the multitude of 
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bands at the Fermi level near the T point. While the 
topologically nontrivial superconductivity in such 2DEGs 
has been studied with relatively simple models [2H GSj, 
these studies rely on the assumption that only one or two 
bands crossing the Fermi level, which, however, seems to 
be at variance with the reported experimental data [26] . 
This is because, unlike in the simplified models, the con¬ 
duction bands of these 2DEGs near the T point cannot 
be attributed to a single set of the transition metal t 2 g 
orbitals. Given that the formation of 2DEG requires a 
confining potential, it is likely that multiple quantum 
well channels arise from each t 2 g orbital [27] , Exactly 
how many bands occur at the T point is often diffi¬ 
cult to predict as the 2DEG confining potential is highly 
non-universal. Since the inclusion of more bands can 
turn topologically non-trivial superconductor into topo¬ 
logically trivial superconductor (and vice versa), this is 
an issue that raises question about the robustness of the 
simple model analysis. 

We show here that it is possible to obtain robust Ma¬ 
jorana zero modes in the perovskite oxide 2DEG by us¬ 
ing its crystalline symmetry and anisotropic dispersions. 
Once the crystalline translational symmetry is assumed, 
weak topological indices are well-defined topological in¬ 
variants [25] [25] that are unaffected by how many bands 
cross the Fermi level near the T point. Physically, in 
the case of 2D superconductors breaking time-reversal- 
symmetry, weak topological indices can tell us whether 
there would be a Majorana zero mode on an edge dis¬ 
location nasi], which can be detected with the STM 
tip. Although a Fermi surface crossing the boundary of 
the first Brillouin zone (BZ), a requirement for any non¬ 
trivial weak indices, has not been observed yet in the ox¬ 
ide 2DEGs, it is possible with currently available experi¬ 
mental techniques to tune the system to satisfy this con¬ 
dition. This is thanks to one universal feature of the (001) 
perovskite oxide 2DEG - the strong anisotropy of the 
band dispersion. This feature, signified by the sharp dis¬ 
tinction between the light mass and heavy mass bands in 
the (001) perovskite oxide 2DEG ARPES (angle-resolved 
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photoemission spectrum) data |13H15I . means that there 
need not be large changes in either the Fermi level or the 
number of electrons per unit cell in tuning the system 
from the band bottom to the Lifshitz phase transition 
point. Moreover, this Lifshitz phase transition would in¬ 
volve only a single heavy mass band as all the light mass 
bands would be at much higher energy. Indeed, this tun- 
ability makes the oxide 2DEG a unique physical system 
to realize the dislocation Majorana zero mode compared 
to the ones discussed previously [821485] as we shall show 
below. 

II. RESULTS 


two bands, one giving rise to the outer Fermi surface 
closer to the van Hove singularity at crystalline momen¬ 
tum points X = (7T, 0) and Y = (0,7r) and the other giv¬ 
ing rise to the inner Fermi surface closer to the T point. 
Microscopically, the hybridization between d xz and d yz 
orbitals is mainly due to the on-site atomic spin-orbit 
coupling K a so = — A c\s z T y c r and the next-nearest 
neighbor hopping K nnn = -t" XU( c r+e x +e y s o'TxC r - 
c r_g sqt x Ct) +h.c.], where s Q are Pauli matrices with 
spin indices and r Q with orbital indices. From these hy¬ 
bridization terms, we obtain the outer Fermi surface dis¬ 
persion of 

£(k) = — (t + t')(cosk x + cosky) — fi (2) 


A. Band structure of the (001) perovskite oxide 
2DEG 

It has been known that electrons close to the Fermi 
level in typical (001) perovskite 2DEGs are mainly from 
t 2 g bands; the bands formed by d xz and d yz orbitals are 
quasi-ID while the one by d xy orbital is quasi-2D. As 
we are mainly interested in weak topological supercon¬ 
ductivity in such systems, we shall focus on the elec¬ 
trons at the Brillouin zone boundaries, namely k x = ir 
or k y = 7r. One important feature of the (001) per¬ 
ovskite 2DEG band structure is that the electrons close 
to the Fermi level with k x = 7r ( k y = n) originate mostly 
from the quasi-ID d yz ( d xz ) orbitals. This is because 
the low-energy physics at the BZ boundary arises out of 
the heavy-mass bands, and the quasi-ID nature of the 
d xz (d yz ) orbital implies the heavy-mass dispersion in the 
y(x)-direction as well as the light-mass dispersion in the 
a;(?/)-direction. This anisotropic dispersion can be cap¬ 
tured by the tight-binding model, 

KO = — ^[ C r+e x ,cc j( T Cr ,a;,CT + C r+e y ,y,<r C r,y,(T + ^- C -] 

T,f7 

~ ^ [ C r+e y ,x,o- Cr ,a!,cr + c £+g x ,y ,cr C r,y.<T + h-C.] ,(1) 

r,cr 

where cj a a creates an electrons at site r with spin po¬ 
larization a =f, and orbital a = x,y (representing d xz 
and d yz orbitals respectively). Because of the quasi-ID 
natures of d xz and d yz orbitals, \t\ \t'\. This simple 

model is sufficient to explain why we can reach the Lif¬ 
shitz transition by lifting the Fermi level only in the order 
of the heavy-mass dispersion bandwidth 4t' and the elec¬ 
tron filling fraction only by ~ yjt' /i; this would involve 
raising the Fermi level by ~ O.leV and adding 0.3 elec¬ 
trons per unit cell when compared to the KTaOa surface 
ARPES data P2| [3B . Such shift in the Fermi level can 
be achieved by both the electrical gating and optically 
induced oxygen vacancies nanaEziEsi. Meanwhile, 
the contribution from the d xy orbital is suppressed as 
it has the light-mass dispersion in both the x- and the 
y-direction. 

We further need to consider the hybridization between 
the d xz and d yz orbitals in order to obtain from them 


— — t') 2 (cos k x — cos ky) 2 + A 2 + (4 1" sin k x sin k y ) 2 . 

It is clear that, for (t — t') A ,t”, the orbital 
hybridization would have little effect near the X/Y 
points except for shifting the Lifshitz transition to /i = 

— ^/4 (t — t') 2 + A 2 . Even with a large A, as shown in 
Fig{l](a) where we used t = 10t' = 0.5eV and A = 0.26eV, 
the former approximating the first principle calculation 
for KTa 03 T3 S 3?). 30] while the latter larger by roughly 
a factor of 2, the distinction between the heavy-mass and 
the light-mass bands remains sharp. Hence the physics 
near the (7r, 0) point would be dominated by the d yz and 
near (0,7r) by the d xz orbital. 

Since the inversion symmetry is obviously broken in the 
surface 2DEG, the spin degeneracy at the Fermi surface 
should be generically split by the non-zero Rashba spin- 
orbit coupling. For our analysis, it will be sufficient to 
consider only the most generic Rashba term, which is 
orbital independent, 


4°lo = « oE 


cl +& jSyT 0 c r - cl + ~ is x T 0 c r + h.c. 


(3) 


Further discussions on the spin-momentum coupling term 
will be presented in Supplementary Material |IVB| 

After taking into account the hybridization as well as 
spin-orbital couplings, the band dispersions are described 
by K = K 0 + k aS o + k nnn + kpg Q . As the chemical 
potential fi moves, there is a Lifshitz transition at which 
the outer Fermi surface crosses the van Hove points at X 
and Y. As we approach the Lifshitz transition, the low- 
energy band structure near the ( 7 r, 0) point, which would 
mainly originate from the d yz orbital, can be given by 
the first-quantized Hamiltonian 

Ho{k x ,k y ) = —2 1' cos k x — 2 1 cos k y — ji 

—2ao [s y sin k x — s x sin k y ]. (4) 

Likewise the band structure in the vicinity of (0,7r), 
which would originate mainly from the d xz orbital can 
be obtained by the 7r/2 rotation of the momentum and 
the spin in Eq. Q. It is also analogous to the Rashba 
wire that the spin degeneracy at (0,7r) remains unbro¬ 
ken, which means that the Fermi surface splitting does 
not lead to two separated Lifshitz transitions. 
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FIG. 1. The band structure of the (001) perovskite oxide 2DEG near the Lifshitz transition, (a) shows the light- and heavy- 
mass band dispersion (in blue and red, respectively), along k y = n with (solid line) and without (dotted line) the orbital 
hybridization, (b) shows the lower band from the d xz /d yz orbitals after the orbital hybridization, with the spin degeneracy 
removed by the Rasliba spin orbit coupling ao = 0.05eV and the perpendicular Zeeman field h z = 0.05 eV. Note that the 
Lifshitz transition point is split, allowing a single hole pocket without spin degeneracy around (77, 7 r). 


To possibly realize isolated Majorana zero modes, the 
final component needed for the band structure is the Zee- 
man field. Near the Lifshitz transition, there will be both 
higher density of oxygen vacancies near the surface as 
well as enhancement of the quasi-ID characteristics of 
the d xz , yz orbitals. Both can give rise to ferromagnetism: 
the former [38J |4Tj because of the oxygen vacancy acting 
as the magnetic impurity [22] while the latter through 
the inter-orbital Hund’s rule coupling [33]. Both of these 
effects should be amplified by the enhanced density of 
states near the van Hove singularity that occurs at the 
Lifshitz transition. We will consider the ferromagnetic 
ordering in the perpendicular direction as was observed 
in the experiment with the density of oxygen vacancy 
l38j . Then, the ferromagnetism-induced Zeeman cou¬ 
pling K z = hz ^2 r cls z ToC r shall split the Lifshitz transi¬ 
tion into two separated ones, as shown in Fig[l](b), giving 
rise to a finite range of /j, for which there is a single hole 
pocket without spin degeneracy around the M = (7 r, 7r) 
point; for this plot we used t" = t' = 0.05eV with a 0 = 
0.05eV and h z = 0.05 eV. 


B. Dislocation Majorana zero mode in proximity 
induced superconductivity 

For the superconducting state, we will first consider 
the case where the pairing is induced through proximity 
to the conventional s-wave superconductor. This will en¬ 
sure the s-wave pairing in the oxide surface 2DEG. We 
also point out that inducing superconductivity through 
proximity effect can have the advantage of achieving su¬ 
perconductivity at higher temperature. To enhance the 
pairing gap magnitude, we would need strong tunnel¬ 
ing between the superconductor and the oxide surface 
2DEG. This can be achieved through using the higher- 


T c two-band superconductors such as FeSe [22H25] : note 
that by symmetry, the single orbital superconductor is 
unlikely to have a strong tunneling to the both the d xz 
and d yz orbitals. Hence our heterostructure will consist 
of the capping two-band s-wave superconductor on the 
(001) surface of SrTiOa or KTaOa as shown in Fig{2](a). 

The combination of the Zeeman field h z and the s- 
wave pairing gap | A s | in the oxide 2DEG near the Lifshitz 
transition can give rise to the non-trivial weak index, 
v = (1, 1 ), i.e. non-trivial ID topological invariants along 
k x ,y = 7r. For instance, the following low-energy effective 
BdG Hamiltonian with k x = tt is exactly equivalent to 
the Rashba-Zeeman wire superconducting state [H HZ]. 

'kisBdG^kx — T^iky) ( 5 ) 

= Hz [2t(l - cos ky) - S/J. + 2 a 0 s x sin/c y ] - s z h z + /i x |A s |, 

where /r a ’s are Pauli matrices acting on the particle- 
hole Nambu space, Sfj, is the deviation of the chemical 
potential from the value at the Lifshitz transition for 
h z = 0, and we use the basis (ckf, C 14 , — c_ k c! k ^). 
It is well known that this ID BdG Hamiltonian is topo¬ 
logically equivalent to the Kitaev chain [48] (class D 
[25] ) when h 2 z > |A S | 2 + (Sfj,) 2 . The C± v point group 
symmetry of the perovskite surface 2DEG dictates that 
T~L S BdG(k x ,ky = 7 r) should be topologically identical to 

disBdGikx = 7T, k y ). 

The orbital hybridization will not affect the topologi¬ 
cal nature of the superconductivity as long as the s-wave 
pairing is intra-orbital. Given that the s-wave pairing 
has no spin dependence, we see that even if we take into 
account the band hybridization and write diikx = 7r, k y ) 
in the band basis, the Zeeman coupling and the s-wave 
pairing terms will remain unchanged, and hence so re¬ 
main the condition for the topologically non-trivial su¬ 
perconductivity; we leave the detailed discussion to Sup¬ 
plementary Material. 
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(d) LDOS on sSC 



FIG. 2. (a) shows the schematic setup of our system, with the oxide 2DEG superconductivity induced through the proximity 
effect and a pair of edge dislocations is on the oxide (STO or KTO) but not on the s-wave superconductor (sSC); an STM tip 
probes sSC. (b) and (c) shows the wave function profile for the Majorana zero modes on both the oxide surface and the s-wave 
superconductor, (d) plots the local density of states on the s-wave superconductor for various different points of the system, 
with the dark blue curve (point a) being taken right at the oxide dislocation position; note that both the sharp peak right 
above the oxide dislocation and the Majorana zero mode being the only subgap mode below the induced oxide pairing gap of 
~1.2 meV. 


Because of the nontrivial weak topological indices v = 
(1,1), unpaired Majorana zero modes occur at disloca¬ 
tions whose Burger’s vector B in units of lattice spacings 
satisfies B ■ v = 1 (mod 2), where mod 2 is from the 
Zi nature of weak topological indices in class D [50]. To 
confirm this, we have performed BdG calculations of the 
lattice models describing the 2DEG in proximity to a 
two-orbital s-wave superconductor (sSC). As shown in 
Fig{2] (b) and (c), we have obtained the Majorana zero 
mode at each dislocation from the numerical exact di- 
agonalization of the real-space BdG Hamiltonian. Our 
calculation was done on a 240 x 240 unit cell with pe¬ 
riodic boundary conditions. Two edge dislocations with 
the Burger’s vector B = ±e x are placed by one half sys¬ 
tem size in the ^-direction, with the links between the 
dislocations shifted as shown in Fig{2] (a) (see Methods 
for details on implementation). This oxide surface with 
the pair of dislocations is coupled by tunneling amplitude 
of ti = 0.05eV to the s-wave superconductor. The sSC 
has the band structure well-matched with that of the ox¬ 
ide surface (see Methods for the band structure details) 
and the pairing gap of |A S | = 0.05eV. Fig [2] (b) and (c) 
shows the probability distribution of the dislocation zero 
energy states, showing sharp peak for both the oxide sur¬ 
face and the sSC, even though the latter does not have 
any dislocation. 

This wave function profile suggests that the STM 
would be a good experimental probe on our disloca¬ 
tion Majorana zero mode 5T|. When the STM tip is 
brought to the sSC as shown in Fig. [ 2 ] (a), the lo¬ 
cal differential conductance Jp(r, uj) is proportional to 
the local density of state (LDOS) of the sSC, p(r,w) = 
J2i[\ u i( r )? S (u ~ Ei) + k(r)| 2 (5(w + Ei)] where the u^Vi 
are the electron and hole components of the i-tli en¬ 


ergy eigenstates, up to replacing the delta function by 
a Lorentzian with the width given by the STM energy 
resolution, which is chosen to be O.lrneV for Figj2] (d). 
We therefore predict that the STM will see a sharp zero 
bias anomaly when it is brought to the point on the sSC 
that is right over the dislocation, the point a of Figj2] 
(d). This anomaly is unambiguously separated from the 
signal of other low lying states, which has a minimum 
energy of the induced oxide bulk pairing gap ~ 1.2meV. 
This is because the Majorana zero mode is the only 
midgap state localized at the dislocation, unlike at the 
Abrikosov vortex where other low energy (~ , where 

A is the pairing gap and Ep the Fermi energy) bound 
states are present. Hence the zero bias anomaly in the 
crystalline dislocation can be regarded as more unam¬ 
biguous signature of the Majorana zero mode than that 
of the Abrikosov vortex [52] , 


C. Dislocation Majorana zero mode in intrinsic 
superconductivity 

We now consider the case of intrinsic superconductiv¬ 
ity in the oxide 2DEG without proximity to conventional 
superconductors. When the oxide 2DEG becomes su¬ 
perconducting at this electron density, there arises possi¬ 
bility of a protected Kramer’s doublet of Majorana zero 
modes at each dislocation when no Zeeman field is ap¬ 
plied. Due to the Rashba spin-orbit coupling, the in¬ 
trinsic superconductivity should generically have on the 
Fermi surface a mixture of the s-wave pairing and the 
p -wave pairing, the latter with momentum-dependent 
Cooper pair spin state [53] . For simplicity, we assume the 
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FIG. 3. (a) shows the energy gap A when we have mixture of the time-reversal invariant p- wave and the s-wave pairing, with 
p = 0 giving us the purely p- wave pairing and p = 1 giving us the purely s-wave pairing. Note that for hz ^ 0, there need not 
be any gap closing in going from the pure p-wave to the pure s-wave pairing, suggesting that the two topological phases for 
the p-wave pairing - one at high hz (in pink) and the other at low hz (in blue) - are identical to those of the s-wave pairing. 
By contrast, for the case hz = 0, where the time-reversal symmetry is preserved, the gap closing around p = 0.7 shows that 
the p-wave pairing (in orange) and the s-wave pairing are topologically distinct, (b) shows the probability distributions for the 
two dislocation zero energy states obtained for hz = 0 and p = 0.4, which have exactly identical profile, (c) shows the energy 
level spacing for the eigenmodes with non-negative energies, again for hz = 0 and p = 0.4. For all positive energy eigenstates, 
separated from the zero energy by a gap of ~12meV, there is double degeneracy due to the time-reversal symmetry. At the 
zero energy, we used the full and dotted lines to indicate the quadruple degeneracy at the zero energy from the occupancies 
and vacancies of the two zero energy states originating from the Kramer’s double of Majorana zero modes at each dislocation. 


following intra-orbital on-site and the nearest-neighbor 
pairing preserve the time-reversal symmetry [54] : 

Hpair ~ |A S | ^ ' ( C r.t 7 b C r ..l. 4“ h.C.) 
r 

~l A tlX! [4+e,( s i^) T o4 +h.e.]. (6) 

r i=x,y 

With this pairing, the BdG Hamiltonian along the k x = n 
cut, 

T-L t BdG{kx = 7r) = /Uz[2£(l - cos k y ) - 5p + 2 a 0 s x sinfc^] 
+/jl x (\A s \ - lAtlsjcSinfcj,), (7) 

will be that of a time-reversal invariant ID topological 
superconductor in class Dill [50], [55] when 5p > 0 and 
|A S | < |A t ||sin/c y | is satisfied at the Fermi surfaces so 
that the gaps at the two Fermi surfaces have the opposite 
signs. In that case, there exists of a branch of helical 
Majorana edge state around k X}V = tt. This means that, 
when we use the argument of the previous section with 
the additional constraint of the time-reversal symmetry, 
there should be a Kramer’s doublet of Majorana zero 
modes at a dislocation with the Burger’s vector of B = x 
or y. 

When the Zeeman field is non-zero, there can be a 
“re-entrant” unpaired Majorana zero modes at the dislo¬ 
cation. This is because Eq. (JtJ) with the addition of the 
Zeeman field 

T~L t BdG{k x = 7r, k y ) = [2f (1 - cos k v ) - Sp + 2 a 0 s x sin k y ] 

Sz^z T Mx(|A s | | A^|s x sin k y ), (8) 

in the h 2 z > (Sp) 2 + |A S | 2 regime is topologically equiv¬ 
alent to the topological phase of Eq. §. Hence, for 


Sp > 0 at hz = 0, if the triplet pairing dominates, 
i.e. |A s | < |A t ||sinfc y | at the Fermi surfaces, there is 
a Kramer’s doublet of Majorana zero modes at each dis¬ 
location for hz = 0 and a single Majorana zero mode for 
^>(<5/r) 2 + |A s | 2 . 

Given that gap closing cannot be avoided at the topo¬ 
logical quantum phase transition, Fig. [3] (a) gives us 
the complete topological phase diagram for the intrin¬ 
sic superconductivity at a fixed Sp = 0.02eV. In this 
plot, the tuning parameter p is introduced to determine 
the relative strength of the p -wave and s-wave pairings, 
i.e. A = pA s + ti(1 — 77) | A t o|(s y sin k x — s x sin k y ) where 
A s = |A s0 |[ti + T2(2/7t)(1 — 77) a,ictan(hz/hzo)] (we have 
set | A s0 1 = 0.04eV, |A t0 | = 0.08eV and hzo = 0.01 eV) 
[56 . For hz 7 ^ 0, as shown in Fig. [3ja), it is always 
possible to adiabatically tune p from 0 to 1 without gap 
closing, while for any value of p , one cannot increase hz 
without closing the bulk gap at some point. We therefore 
conclude that the non-trivial high (trivial low) hz phase 
at p = 1 is topologically equivalent to that of p = 0 , the 
purely s-wave pairing case we considered in the previous 
subsection. However, if we restrict ourselves to the case 
with time-reversal symmetry, i.e. hz = 0, Fig. [U ( a ) 
shows there is a gap closing around p = 0.7, consistent 
with the fact that the time-reversal invariant p-wave pair¬ 
ing (p = 0 ) is topologically distinct from the pure s-wave 
pairing (p = 1). At p = 0.4 and hz = 0, the two zero en¬ 
ergy states localized at the two edge dislocations with ex¬ 
actly identical probability distribution, as shown in Fig. 
[3])b), indicates the existence of the Kramer’s doublet of 
Majorana zero modes at each dislocation. This confirms 
that the existence of the Kramer’s doublet of Majorana 
zero modes at each dislocation characterizes this time- 
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reversal invariant topologically non-trivial phase. This 
means that with an STM with an s-wave superconduct¬ 
ing tip over this dislocation, we should be able to observe 
time-reversal anomaly m- As in the case of the prox¬ 
imity induced superconductivity, these dislocation Ma- 
jorana zero modes, as shown by Fig. [3j^c) , are the only 
subgap modes of the system. 


III. DISCUSSION 

We have shown in this paper how isolated dislocation 
Majorana zero mode can arise from both the proxim¬ 
ity induced and intrinsic superconductivity in the oxide 
2DEG. Its existence can be considered the most perti¬ 
nent criterion for the topologically non-trivial supercon¬ 
ductivity in the oxide 2DEG, and it can be experimen¬ 
tally detected through STM. The crucial requirement for 
achieving such superconductivity is that the oxide 2DEG 
needs to be close to the Lifshitz transition. 

The key difference between the proximity-induced and 
the intrinsic oxide 2DEG superconductivity is that the 
Zeeman held is a necessary condition for the non-trivial 
topology in the former but not for the latter. The physi¬ 
cal consequence is that for the intrinsic superconductivity 
in the absence of the Zeeman held, the dislocation can 
host a Kramer’s doublet of Majorana zero modes; this is 
not possible if the superconductivity is induced through 
proximity to an s-wave superconductor. By contrast, in 
the presence of nonzero Zeeman held, the only possible 
protected midgap state on a dislocation is a single Majo¬ 
rana zero mode regardless of the origin of superconduc¬ 
tivity. 

While our intrinsic superconductivity with the non¬ 
trivial weak index at the zero Zeeman held has the es¬ 
sentially same pairing symmetry as the topological su¬ 
perconductivity investigated by Scheurer and Schmalian 
[25] , these states are topologically distinct. From Eq.([7|, 
our non-trivial phase requires Sy > 0 while that of 
Scheurer and Schmalian requires Sy, < 0 with the Fermi 
surfaces enclosing the T point, and with this pairing sym¬ 
metry the gap closing around <5/x = 0 cannot be avoided. 
This reflects the fact that, with the reflection symmetry, 
the topological invariant of the Dill class in 2D can be 
Z rather than Z2 [551E3] ■ The existence (absence) of the 
dislocation Majorana doublet for 5y > 0 (Sy < 0) can 
be regarded as a physical manifestation of this topolog¬ 
ical distinction. We leave to future work what type of 
interaction would favor this pairing symmetry near the 
Lifshitz transition. 

Lastly, we want to point out that it is generically easy 
to change the topology of the superconducting state of 
the (001) perovskite oxide 2DEGs. This is because the 
universal anisotropic band structure makes it easy to ac¬ 
cess the van Hove singularity through gating and opti¬ 
cally inducing oxygen vacancies. While there have been 
previous works on the physical realization of the the dis¬ 
location Majorana zero mode [52H55] , they have not pro¬ 


vided easy means to alter the weak indices of the su¬ 
perconducting states. Therefore we conclude that not 
only is the dislocation Majorana zero mode the most ro¬ 
bust topological feature of the oxide 2DEG superconduc¬ 
tor but also that the oxide 2DEG superconductor is the 
particularly suitable system for realizing the dislocation 
Majorana zero mode. 


IV. METHODS 

A. Weak indices and dislocation Majorana zero 
modes 

It is possible in a 2D superconductor on a square lattice 
to consider the ID topological invariants defined along 
k x ,y = 7r, which are known as the weak indices |60ft62) . 
In general, the weak index z/j can be defined for each 
time-reversal invariant momentum G,/2 (which makes 
G, a reciprocal lattice vector) as a topological invariant 
of the manifold perpendicular to G, but contains Gj/2, 
and hence the weak indices can be written as a single 
vector v = z^Gj, where G, is the unit vector parallel 
to G, . The C\ v symmetry of of our 2DEG means that 
its u will have only a single independent component v 
and therefore can be written as v = i/(l, 1). The weak 
indices is clearly topologically protected when the sys¬ 
tem has crystalline symmetry, the topological crystalline 
insulators[63] being one class of examples. In this paper, 
we will focus on its manifestation through the Majorana 
zero mode localized at its crystalline topological defect - 
the edge dislocation [5UH55] , 

We first note that the non-trivial weak indices in a 
superconductor imply the existence of a branch of Ma¬ 
jorana edge modes around k e d ge = n. Since restrict¬ 
ing ourselves to the k x = n manifold means converting 
the 2D mean-field Hamiltonian HsdG{k x , k y ) into the ID 
Hamiltonian HBdc(k x = 7r, fc y ), the non-trivial weak in¬ 
dex means that, for the simplest case of the class D, where 
the time-reversal symmetry is broken, a single protected 
Majorana zero mode exists at k x = 7r for the edge run¬ 
ning in the x-direction. This is possible only if there is 
a branch of chiral Majorana edge state centered around 
k x = 7 r. Note that the existence of this branch of the 
edge state is determined by the projection of v to the 
time-reversal invariant momentum (zr, 0). 

A single Majorana zero mode exists at the edge dislo¬ 
cation when there is a chiral Majorana edge state cen¬ 
tered around k e d ge = 7r - To see how this arises, note that 
the dislocation can be created by severing all links, both 
through hopping and interaction, between two halves 
(y < 0 and y > 0) and then non-trivially re-connect 
the two halves to introduce the edge dislocation, with 
the x < 0 part glued back according to the original links 
but the x > 0 part has all the links altered by trans¬ 
lating the sites of the y > 0 half by a lattice constant 
along the x-direction, which sets the Burger’s vector of 
this dislocation to be B = x [25] HD]. Now when this 
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system was cut, there would have been Majorana edge 
states along the x-direction for both y < 0 and y > 0 
with the opposite chirality. Hence when the system is 
glued back along the original links, the tunneling between 
the two edges would lead to the backscattering that gaps 
out these edge modes, with the mass gap being propor¬ 
tional to the tunneling amplitude. However, when the 
dislocation described above is introduced, there will be 
a qualitative effect on the tunneling between the k x = ir 
edge state. This is because the k x = 7r edge mode wave 
function reverses its sign when we translate by one lattice 
site along the x-direction, the relative sign of the k x = tt 
edge modes for y < 0 and y > 0 edges will change its 
sign at the dislocation. That means that if a dislocation 
is introduced when we glue back with only infinitesimally 
weak coupling across y = 0, the effective low energy ac¬ 
tion along y = 0 for the k x = 7r edge modes would be 


Seff — 



id t + id x im 0 sgn(x) 
—imosgn(x) idt. — id x 




(9) 


where the upper and the lower component correspond 
to the upper and the lower edge and m o is proportional 
to the tunneling amplitude for the k x = tt modes; this 
action is well-known for having a single Majorana zero 
mode at our dislocation x = 0: 


r 0 = 


y/rno 

2 


exp(-m 0 |x|) 


1 

1 


( 10 ) 


By contrast, the existence of the k x = 0 branch is irrele¬ 
vant as its tunneling amplitude does not change sign at 
the dislocation. Since the Majorana zero mode is pro¬ 
tected as long as it remains separated from other Majo¬ 
rana zero mode, the Majorana zero mode that arose at 
the infinitesimal coupling across the y = 0 cut will persist 
when the coupling across y = 0 is increased to the bulk 
values. In general, the condition for the existence of the 
protected Majorana zero mode is u ■ B = 1 (mod 2). 

We can similarly show the existence of a Kramer’s dou¬ 
blet of Majorana zero modes at the edge dislocation when 
there is a helical Majorana edge state centered around 
kedge = 7r ■ The key point here is that dislocation involves 
no time-reversal symmetry breaking and therefore, in the 
‘cut and paste’ picture, the Kramer’s doublet needs to 
be maintained even with the inter-edge backscattering. 
Therefore, when is y is the intra-edge time-reversal op¬ 
eration, the effective low energy action for the k x = tt 
helical edge mode would be 


S e ff — / dx 




idt + is z d x 
—i(s z m + s x m') sgn(x) 


i{s z m + s x m')sgn(x ) 
id t - is z d x 


*■( 11 ) 


This action gives us two Majorana zero modes, 


r+ = 


y/rfio 


exp(— m\x\) 


cos to | a; | 
sin m! I x I 


cos to | a; | 
sin to'I a; I 


( 12 ) 


T_ = exp(—m|x|) 


— smm \x\ 
cos m'| a; | 

— sin to 7 | a: | 
cos to' lx I 


(13) 


which form a Kramer’s doublet, i.e. (is^jTj- = Trip. 


B. Real space Hamiltonian with dislocation 

We need to have the real-space BdG Hamiltonian in 
order to obtain the dislocation Majorana zero mode 
through exact diagonalization. We first note that the 
terms in our real-space Hamiltonian could be divided into 
three groups, the first being the onsite term, 

H-onsite ^ ) ( M T a a S[)TQC r a r j 

r 

A z TyC r H - h.c.) + H ons n e 

—pair > (14) 

r 

consisting of the chemical potential, the Zeeman energy, 
the atomic spin-orbit coupling and the onsite pairing, the 
second being the nearest neighbor terms 

^(r), 

r i=x,y 

Hg.(r)=-t ^ ^(Cr+ei.a^OCr.a + h.C.) 

a—x,y 

^ ^ ^ (1 (c r +ei,a^0^ r > a 

a=x,y 

+ a 5Z[^+ei,a(^ 5 i) C r,a + t.C.] 

a 

H - H nn —p a i r (r, (15) 

which includes the spin-conserving intra-orbital near¬ 
est neighbor hopping, the Rashba spin-orbit coupling, 
and the nearest-neighbor time-reversal invariant pairing. 
Lastly, we have the next-nearest neighbor hopping, 

Hjinn — EE Ke x +se y ( r )? 

r s=± 

(r) t" Y: E + h.c.). (16) 

r s=± 

which gives us the spin-independent component of the 
orbital hybridization. For the proximity-induced super¬ 
conductivity, we set the nearest neighbor pairing to be 
zero, i.e. H nn _ pa i r (r,ei) = 0, and set the onsite pairing 
to be originated entirely from a two-band s-wave super¬ 
conductor: 

Honsite—pair ~ L ^ ( (CpS'QTQ./ r -f h.C.) -f H s , (17) 



where H s is given by 

H s = t s 'y ^ ^a,i(/ r _|_£ ija 5o/r,a “1“ h.C.) 

a=x,y 

-t's E (! - <M(/E i ,a 1 2 3 4 S 6 7 o/r,a + h.C.) 
a=x,y 

- Ms E /rWo/r + |A S |E(/ r Vo/h + h - C -)( 18 ) 

r r 

(in Fig|2j we have set t s = t, t' s = t' and p s = —0.7eV 
for the maximal proximity effect), while for the intrinsic 
superconductivity calculation shown in Fig|3j we set the 
pairing terms to be 

H onsite—pair =77|A s0 |[1 - i{ 2/7t)(1 - rj) arctan(/i z /ft zo )] 
x E( C r,t T ° C r 4 + h - C -)’ 

r 

H„n-pair{r,ei) = - (1 - 77)|A t0 1[4 + g. (s i s a ,)r 0 c); + h.c.], 

(19) 

where we have set |A s0 | =0.04eV, |A t0 | =0.08eV, 
h-zo =0.01eV with 0 < rj < 1 determining the relative 
contribution of the p- and s-wave pairings. For the in¬ 
trinsic superconductivity under a finite hz, the phase 
difference between the s-wave and the p -wave pairing 
were added so that the pairing terms would break time- 
reversal symmetry. 

In the real space, the dislocation point serves as a 
starting point for a branch cut along which the nearest- 
neighbor Hamiltonian H ^ (r) is applied on a next-nearest 
neighbor link. Our square lattice N x = 240 by N y = 240 
latitce has a periodic boundary condition to both direc¬ 
tion. In order to have a dislocations at {N x /4, N y / 2 + 1) 
with the Burger’s vector B = +e x and another dislo¬ 
cation at (3N x /4, Ny/2) with B = — e x , we apply He y 
on the links connecting (n, N y / 2) and (n + 1, N y /2 + 1) 
for N x / 4 < n < 3N x /4, Ke l: +e y on the link con¬ 
necting (n,N y / 2) and (n + 2,N y /2 + 1) for N x / 4 < 
n < 3N x /4 — 1, and K^-ey on the link connecting 
(n, Ny/2 ) and (n, N y /2 + 1); meanwhile between the two 
nearest neighbor pairs (N x /4,N y /2) and (N x /4, N y /2 + 
1), (3N x /4, Ny/2) and (3N x /4,N y /2 + 1) and also be¬ 
tween the two next-nearest neighbor pairs [N x /4, N y /2 + 


1) and {N x /4 + l,N y /2), {3N x /4 - l,N y /2 + 1) and 
(3N x /4, Ny/2), all hoppings and pairings are set to zero. 
Note that for the case of proximity-induced superconduc¬ 
tivity, the s-wave superconductor remains completely free 
of crystalline defects. 

While we set some of the parameters to be rather large 
for the sake of convenience in the numerical calculation, 
such choice does not affect the topological properties of 
the system. For instance, A = 0.26eV is about factor 
of 2 larger than the estimated value for the tantalum 
atom, while a 0 = 0.05eV is several times larger than the 
estimated value from the magnetoconductivity measure¬ 
ment [23]. These choices are intended to increase the 
bulk energy gap so that our lattice size is sufficient to 
see a localized dislocation zero mode. This increase in 
the bulk energy gap occurs away from the BZ bound¬ 
ary k XyV = 7r, e.g. the larger a o increases the energy 
gap along k x = ±k y , while the larger A lifts the higher 
d X z/yz band away from the Fermi level. Such changes 
do not affect weak indices, which are the ID topologi¬ 
cal invariant along k x , y = n. Concerning real materials 
and experiments, as long as the system is in the topo¬ 
logical regime and the induced bulk gap is large enough, 
i.e. much larger than the STM resolution, the disloca¬ 
tion Majorana zero mode and the zero energy anomaly 
can be detected clearly as shown in Fig{2jc). 

Lastly, we point out that with our p -wave pairing in 
Fig. i (a) for hz = 0 allows for a finite range of r) for 
which there are nodal quasiparticles. While it is possible 
in principle to come up with a p-wave pairing for which 
the energy gap closes for only a single value of r), such 
p -wave pairing should have constant magnitude over the 
entire Fermi surface, which in general is not possible with 
our nearest-neighbor pairing. 
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SUPPLEMENTARY MATERIAL 


A. Rashba spin-orbit coupling near van Hove singularity 


The C± v point group symmetry of the (001) perovskite surface 2DEG allows for the spin-momentum coupling in 
the d xz , yz orbitals from the following nearest-neighbor hopping terms 12511661 : 


k-RSO — 'y ' ^r+zjSyToC, 




— c. 


r+By 


is x TQC T + h.c. 


+ ai 


E[« 


'i'Sx'l’xCr C r j.s iSyT x C r h.C. 


+ OL 3 


£[■ 


'r+ea; 


C „_LA ^SyT Z C r + C. 


r+e, 


iSxTzC r + h.C. 


(SI) 


This raises the question how these terms arise microscopically. Also we can ask whether there is any qualitative 
effects due to the orbital dependent terms 0 : 1 , 3 , i.e. whether there is any sin fc^-dependent spin-momentum locking 
of magnitude in the order of 00 , 1,3 along k x>y = tt. 

The starting point for the microscopic physics is that due to the breaking of the inversion symmetry with respect 
to the xy- plane, there can be spin-independent hopping between d xz , yz orbitals and other d-orbitals: 


K inv = ~ 7l y! Xl 1 - 8a,a') <4+e a ,a' S 0Cr - c[J + g a S 0 C r,a' + h - C - “ 72 X X ^ a - a ' C i!-+e a ,a' S QCr - 4+e a S 0 C r , a , + h.C. 


- 73 X) [ c i+e x ,x s oc , r - cJ^SQCr,* - cl + ~ y y s 0 c' r + g^SoCr,* + h.c. 


(S 2 ) 


where c!,c,c! are the annihilation operators for d xy ,d z 2 ,d x 2_ y 2 orbitals respectively. Note that, in the bulk, any 
spin-conserving hybridization of this type would not be allowed due the fact that these orbitals are even under the 
inversion with respect to the xy -plane while the d XZtVZ orbitals are odd. But even in the bulk, this hybridization is 
allowed through the transition metal spin-orbit coupling: 


K'aSO ~ ~ 


[4,x(-* S *) C r + C i,y(™ y )Cr + h-C. 

r 

V^3A ^ ^ j^r,x( iSy^c r T c][ y (zs a .)c r -t- h.c. A ^ ~ [c\ x {is y )c T -\- cj y {is x ^c r T h.c. 


(S3) 


The hybridization of the d XZtyz orbitals with the other d-orbitals through K lnv + K' aSO effectively generates the 
Rashba-Dresselhaus terms of Eq.|Sl|) in the subspace of d xz ^ yz orbitals [251 ESI [67]- One relatively simple way to 


obtain these Rashba-Dresselhaus terms is through the second order degenerate perturbation theory 


,v\fcRSo\d',cr') = X 
b^x,y 


(a,a\K inv \b, a")(b , er"|A^ so |a / , a') + (a, a\K' aSO \b, cr")(b, a"\K mv \a', a'} 


£o - 


(S4) 


the point here being that we take a, a' to be the band (with spin degeneracy) formed by d XZtVZ - 

One question that arises here is whether the orbital b in Eq.(S4) should include the e g orbitals as well. Much of the 


analysis for the Rashba-Dresselhaus term near the T point excludes the e g contributions to Eq. (S4) 23 [Ml ESI EZ] , 


which would be justified in the limit where the e g orbital energies are much higher than that of the d xz , yz orbitals. 
Applicability of this limit is qualitatively important for the analysis near the X point, for excluding the e g contribution 
effectively gives us 


k-RSO = 


4yiA 


s y sin k x 


xy,yz 


(S5) 
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■, xy ,yz is the energy splitting between the d xy and d yz orbitals at the X point, for the d. 


where St;, 

straightforward to obtain analogous result for the d xz 
justifiable to take the spin-momentum coupling to be K^ so of Eq.([3]). However, when the inversion symmetry is 
broken, the energy splitting between the two e g orbitals actually may be in the same order of magnitude as the crystal 


yz orbital (it is 
orbital near the Y point), in which case it would not be 

(o) 


field splitting between the t2 Q and e g orbitals 
from the d z 2 orbital) to Eq.(S4) would give us 


In such case, the contribution of the lower e g band (mostly 


K-rso 


4yi A 

xy,yz 


s v sin k x + 


4a/37 2 A . 

—- s x sin k y . 

,yz 


(S6) 


Hence we are justified qualitatively in setting K^so °f Eq.pl to be the spin-momentum coupling in the d xz yz orbitals 

m- 


B. Orbital mixing effect near van Hove singularity 


In this Section, we will show that even with a strong orbital hybridization, Eq.([6]) still gives us the effective low- 
energy Hamiltonian at k x = tv. 

As a starting point, let us consider the first-quantized Hamiltonian in the normal state for the d xz /d yz orbitals, 


n 


( 0 ) 

normal 


= — (t + t')( cos k x + cos k v ) — n — r z {t — t')( cos k x — cos k y ) + r y s z A + T x At" sin k x sin k y 


(S7) 


(where cq’s are the orbital Pauli matrices, with the d xz / yz being the eigenstate of cr 3 with the eigenvalue of ±1), where 
we have left out the Zeeman and the Rashba terms. When A is in the same order of magnitude as t — t 1 , as is known for 
KTa0 3 , there would be a considerable orbital hybridization even when k x ^ y = 7 r. Along k x = 7 r, we can relate the band 
basis and orbital basis through the transformation U = exp(z<Tis z /3fc B /2) where tan /3k y = —A/[(t — t')(l + cosk y )], 
giving us 


UH 


( 0 ) 

normal 


k x — 7 T 


= (t +t')( 1 — cos k y ) — fi + t z ^J\t — f') 2 (l + cos k y ) 2 + A 2 
« [t +1')( 1 — cos k y ) — 1 + r 


y/4 (t ~ t') 2 + A 2 


(S8) 


The low-energy Hamiltonian will then be given by the projection to the lower band, cr 3 = —1. 

We then need to consider how the other terms, which are much smaller in the magnitude, transforms under U. 
Since it is obvious that the Zeeman term —s z hz and the intra-orbital s-wave pairing T!|A s | remains invariant, we can 
mainly focus on the spin-momentum coupled nearest neighbor hopping. From Eq.(Sl|, we can see that these terms 
in the first-quantized form comes out to be 


ICso = — 2ao(Sa; sin k y — s y sin k x ) + 2ot\T x (s x sin k x — s y sin k y ) + 2 a 3 r z (s x . sinfc y + s y sin k x ). 
In the band basis, these terms comes out to be 

U /Cso| fcx=7r U' =- 2(a 0 cos /3 ky - a 3 T z )s x sin k y - 2T x (a 0 sm/3 ky + a^SySmk. 

2 (t - t') 


(S9) 


' - 2 


a 0 - 


Q 3 t z 


s x sin k y — 2 t x 


A 


a 0 


V4 (t - t'Y + A 2 

We can now see that the t z = — 1 projection of the full BdG Hamiltonian is 

(t-t'f ^ 


y/&(t - ?y- + A 2 


a 1 


Sy sin ky. 


(S10) 


PU U B dG\ k -~ U ] P=v 


{t + 1') 


y/(t — t') 2 + A 2 /4 


(1 — cos ky) — 5[i — 2 < ao 


C t-t') 


V(t-t') 2 + A 2 /4 


+ a 3 > s x sin k y 


- s z h z + //xIAf, 


(Sll) 


where P = (1 — t z )/2 is the projection operator to the lower band; this is clearly in the same form as Eq. §. 






























